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Transport properties in the presence of magnetic fields are numerically studied for the spin- 
1/2 Heisenberg XXZ chain. The breakdown of the spin-reversal symmetry due to the magnetic 
field induces the magnetothermal effect. In analogy with the thermoelectric effect in electron 
systems, the thermomagnetic power (magnetic Seebeck coefficient) is provided, and is numeri- 
cally evaluated by the exact diagonalization for wide ranges of temperatures and various mag- 
netic fields. For the antiferromagnetic regime, we find the magnetic Seebeck coefficient changes 
sign at certain temperatures, which is interpreted as an effect of strong correlations. We also 
compute the thermomagnetic figure of merit determining the efficiency of the thermomagnetic 
devices for cooling or power generation. 
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1. Introduction 

During the last two decades, strongly correlated sys- 
tems with reduced dimensions have been extensively 
studied from both theoretical and experimental perspec- 
tives, due to their unique static or dynamical proper- 
ties. Focusing our attention on the transport properties 
of magnetic materials described by one- or quasi one- 
dimcnsional spin-1/2 magnets, we find their unconven- 
tional features induced by magnetic excitations. For in- 
stance, an unusually large spin diffusion constant has 
been observed in an NMR experiment for S^CuOa 1 
which is well characterized by the spin-1/2 Heisenberg 
chain. In addition, anomalously enhanced thermal con- 
ductivities were also measured in experiments for Heisen- 
berg chain compounds SrCuC>2 and Sr 2 Cu03, 2 ' 3 indicat- 
ing a quasi-ballistic heat transport carried by the spinon 
or magnon with a mean-free path much larger than the 
correlation lengths. 4 Theoretically, these quasi-ballistic 
transport might be interpreted as a residual effect of in- 
tegrability. In fact the energy current of a certain class 
of integrable systems including the spin-1/2 Heisenberg 
chain is written as a constant of motion, 5,6 and there- 
fore the heat transport exhibits purely ballistic behav- 

5-8 



\ ior. 8 On the other hand the spin current is generally 



not conserved even in integrable systems, nevertheless it 
might have a finite overlap with the conserved quantities 
underlying integrability. Consequently the spin current- 
current correlation does not decay to zero for long times, 
and then the spin transport is also considered to be bal- 
listic. 9-12 Although real materials such as Sr 2 Cu03 have 
non-integrable terms as a small perturbation, the sys- 
tem still exhibits the (quasi)-ballistic transport proper- 
ties as already mentioned above. The integrability, how- 
ever, is not a necessary condition to guarantee the exis- 
tence of such ballistic transport properties. In fact, the 
quasi-ballistic heat transport was also observed in spin- 
ladder compounds (Sr,Ca,La)i4Cu2404i, 13-16 which can 
no longer be explained by the integrability. The arising 
question of what kind of non-integrable systems do or 
do not keep the (quasi)-ballistic features still remains an 



intriguing open problem. 11, 12, 17-27 

Another crucial problem to be considered in the trans- 
port properties is effects of external fields. By analogy 
with the thermoelectric effect in electron systems, one 
might expect the existence of the magnetothermal ef- 
fect in spin systems; the temperature gradient along a 
sample causes the magnetic field gradient. For the one- 
dimensional spin systems without external fields, how- 
ever, there is no magnetothermal effect since the system 
exhibits the spin-reversal symmetry. In the presence of 
finite magnetic fields, the situation drastically changes; 
the magnetothermal effect indeed arises due to vanishing 
of the spin-reversal symmetry. 8,28,29 Thus for the com- 
plete understanding of the transport properties under 
finite external fields, we must take into account the mag- 
netothermal effect correctly. Quite recently, by utilizing 
the phenomenological approach, the magnetic Seebeck 
coefficient has been analytically calculated in the mass- 
less regimes of the spin-1/2 Heisenberg XXZ chain. 30 
Most significant, in that work, was the prediction that 
the magnetic Seebeck coefficient changes sign at certain 
temperature and for certain interaction strengths. Un- 
fortunately, the result is limited to the low-temperature 
region, due to the lack of conclusive results for the spin 
conductivity of the Heisenberg chain. 

Motivated by this, in this paper, we will discuss the 
magnetothermal effect for the spin-1/2 Heisenberg XXZ 
chain beyond the limitation of the analytical approach. 
Namely using the exact diagonalization method, we eval- 
uate the magnetic Seebeck coefficient for wide ranges 
of temperatures and various interaction strengths. As 
a consequence, we find the magnetic Seebeck coefficient 
changes sign for interaction strengths in the antiferro- 
magnetic regime, when the magnitude of the magnetic 
field satisfies a certain condition. Moreover we evaluate 
the thermomagnetic figure of merit measuring the effi- 
ciency of the thermomagnetic devices. 

The layout of this paper is as follows. In the subse- 
quent section, we briefly present a general formulation of 
the transport properties in spin systems. The transport 
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coefficients and the thcrmomagnetic power are described 
within linear response theory. In section 3, we provide 
the spin and thermal currents for the spin- 1/2 Hciscn- 
berg XXZ chain, and express the transport coefficients 
in terms of the correlation functions among the current 
operators. Numerical results by the exact diagonaliza- 
tion up to 18 sites are presented in section 4. Section 5 
is devoted to the summary and discussions. In appendix, 
we shortly give exact results for the free fermion (XY) 
model. 

2. Transport Coefficients in Spin Systems 

For later convenience, here we briefly provide a general 
formulation of the transport coefficients and the magne- 
tothermal effect in spin systems. 

Let us consider a system with two currents-the spin 
and heat currents-which flow as a result of forces. The 
transport coefficients relate the currents to the driving 
forces, i.e. the potential and temperature gradients. Phe- 
nomenologically these relations may be written in the 

form 28,29 



Js 



Lss L S q 
Lqs Lqq 



-V0 S \ 
-VT/TJ 



(2.1) 



where Js, Jq are the spin and heat currents, respectively, 
W<t>s the potential gradient (typically the magnetic field 
gradient V0 S = — V/i) and VT the temperature gradient. 
The Kubo formula explicitly gives the coefficients Lij 
({i, j} € {Q, S}) in terms of the correlation functions of 
the current operators: 

r/3 



lim Re 

u^O 



/ dtc-' 1 ^-' 1 ^ &\{Ji{-t-i\)Jj), 
JO JO 



(2.2) 

where G {Q,S}, (3 is the reciprocal temperature 

P = 1/T and (• • •) denotes the thermal expectation value 
per site. Using this relation with eq. (2.1), one obtains 
the transport coefficients for the spin system. The spin 
conductivity a is measured under the condition of no 
temperature gradient VT = 0: 



Js = er(-V0 s ), cr = L ss . 



(2.3) 



On the other hand the thermal conductivity k is defined 
when there is no spin current Js = 0: 



J Q = k(-VT), 



K = 



L QQ 



L QS 

Lss 



Finally as another important quantity which deter- 
mines the efficiency of thermomagnetic devices for cool- 
ing or power generation, we define the thermomagnetic 
figure of merit 



ZT - 



T 2 
L QS 



~ L QS 



(2.6) 



3. Current Operators in the Heisenberg chain 

Here we apply the formulae given in the preceding sec- 
tion to the spin-1/2 Heisenberg XXZ chain. The Hamil- 
tonian is defined as 



L ^ L 

7~C = ^ h kk+1 - - o 



k=l 



k • 



lfcfc+1 



J 



°k CT fe+l + "k+l^fc 



A 



(3.1) 



where 



(al ± ia^)/2 and er£ (a e {x,y, z}) arc 
the Pauli matrices associated with the fcth site of the 
chain. From now on, we assume L is even, and impose 
the periodic boundary condition (crj = cr£ +1 ). The cou- 
pling constants J and A together with the magnetic 
field h determine the ground state properties of the sys- 
tem 31 (see Fig. 1). Since the energy spectrum is invariant 
under both the transformations (h,A4 = J^k^k/ty <-> 
(—h,—A4) and (J, A) <-» (—J, —A), we assume h > 
and J > 0. 











massive 








massless / 









(2.4) 

Fig. 1. Phase diagram of the ground state for h > and J > 0. 



Note that here we have imposed the Onsager relation 
Lsq = Lqs- 

On the analogy of the thermoelectric power (Seebeck 
coefficient) for electron systems, we define the thermo- 
magnetic power (we refer to as the "magnetic" Seebeck 
coefficient) 5, which should be measured under the con- 
dition Js = 0: 

V(j> s _ 1 Lqs 



S = — 



(2.5) 



VT T Lss 

The (magnetic) Seebeck coefficient is a crucial criterion 
to clarify the types of carriers. Namely when the sign of S 
is positive (negative), the spin and heat are dominantly 
carried by the carriers with up-spin (down-spin). 



To evaluate eq. (2.2), we define the spin and energy 
currents from the continuity equations for the local spin 
S% = a^/ 2 and energy h kk+1 : 5 



-divj 



ifcfc+l 



(3.2) 



where divj k := j k+1 — j k . These relations give the spin 
and energy current, Js = Efcife = Efc 'A h k-ik, and 
Je = Efeife = Efci^fe-ifc./ifefe+i]) respectively. Their 
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explicit forms read 

L 



fe=i 



k=i l 

- H°k-1 + <t| +2 )(c7+ - cr+ +1 cr-)|. (3. 



Note that the heat current Jq should be defined as 



Jq = iy^[fefc-i,fc - hS k _ 1 ,h k ,k+i - hSl] = Je 



hj s . 
(3.4) 



In general, the exact evaluation of the dynamical cor- 
relation functions such as eq. (2.2) is highly challenging 
problem even in exactly solvable models. However, the 
energy current of a certain class of solvable models in- 
cluding the present system is written as one of the non- 
trivial conserved quantities underlying integrability, 6 i.e. 
[H, Je] = 0. This important fact directly leads to the 
diverging thermal conductivity at zero frequency uj = 0: 



k = itDq5(uj), 



(3.5) 



where the weight of the delta function Dq is referred to 
as the thermal Drude weight 7,8 ' 28-30 



D Q = p 2 {ji)-0- 



(JeJs) 2 
D S ' 



(3.6) 



Here Ds is the Drude weight for the spin transport, 

a = ttD s S(uj) +(T rcg . (3.7) 



For zero magnetic field h = 0, the magnetothermal ef- 
fect, which is related to the off-diagonal dynamical cor- 
relation Lq S (~ Lsq), is always zero (i.e. (JeJs) = 0) 
because the system exhibits the spin-reversal symmetry. 
In this case the thermal Drude weight is simply written 
as Dq = f3 2 (J^), and can be exactly evaluated by the 
Bcthe ansatz technique. 7,8 

In contrast, for finite magnetic fields h > 0, the mag- 
netothermal effect arises due to vanishing of the spin- 
reversal symmetry, (JeJs) > 0. In addition to this, the 
spin conductivity also diverges (i.e. Ds > 0) as far as 
h > 0, which is proven by using the Mazur inequality 6 



(3-S 



Hence the thermal Drude weight (3.6) and the magnetic 
Seebeck coefficient (cf. (2.5)) 



„ 1 f (JeJs) 1 , 
S -f\-^ s -T- h 



(3.9) 



are both finite at finite temperatures and magnetic fields. 
Insertion of eqs. (3.5), (3.7) into eq. (2.6) yields 



zt = ^Lt. 



D 



Q 



(3.10) 



To analyze the magnetothermal effect for the XXZ 
chain (3.1), we must explicitly determine the spin Drude 
weight as well as the correlations (j£) and (JeJs)- In 



fact, both the quantities (JeJs) an d (Je) have already 
evaluated by the Bcthc ansatz in rcf. 30. On the other 
hand, for the spin Drude weight, we may derive Ds by 
considering the finite size corrections of the "string" solu- 
tions to the Bethe ansatz equation. 10,32 In general, how- 
ever, the validity of the results obtained by applying the 
finite size correction to the string solutions arc highly 
questionable under the circumstance that the breakdown 
of the string hypothesis is reported as in ref. 33-35. In- 
deed the resultant Drude weight Ds derived in ref. 10 
for the case of zero magnetic field is inconsistent with 
that obtained by numerical analysis 11 or field theoretical 
arguments. 12 Consequently the exact analysis of Ds at 
arbitrary temperatures still leaves an open problem. 

Alternatively, in the next section, we will calculate the 
magnetic Seebeck coefficient S and the figure of merit 
ZT by the exact diagonalization for finite chains. 

Finally we comment on the transport properties of the 
spinlcss fermion system given by performing the Jordan- 
Wigncr transformation on the XXZ chain (3.1): 

Hc = Xj{*( c fc Cfe+i + c I+i c fc) + v ( n k - \) K+i - \) J 



L 

fe=i 



(3.11) 



where cl (cfe) is the fermionic creation (annihilation) op- 
erator on the site k and n k — c\ck - Note that we have also 
transformed J — > t, h — > \i and A — > V/(2t). The cor- 
responding charge and energy currents are respectively 
given by 

L 

Jc = ety^icj^Cfc+i +h.c), 

L 



fe=i 



J E = t t(i4_ 1 Cfe+i + h.c.) 



+ V(n k+1 + n k+2 - l)(i4c*+i + h.c.) j, (3-12) 

where e is the charge of the particle. Accordingly the 
heat current is defined as Jq = Je — (fi/e)J c (cf. (3.4)). 
In the thermodynamic limit L — » oo where the difference 
of boundary conditions between the spin system and the 
corresponding fermion system vanishes, the charge (ther- 
mal) conductivity a c (k c ), the Seebeck coefficient S c and 
the thermoelectric figure of merit Z C T are respectively 
related to those of the spin system: 

Z C T = ZT. (3.13) 



e 2 a, 



4. Numerical Analysis for magnetic Seebeck Co- 
efficients and Figures of Merit 

In this section we evaluate the magnetic Seebeck coef- 
ficient S and the thermomagnctic figure of merit ZT by 
the exact diagonalization up to 18 sites. 

Before we analyze them in detail, let us shortly review 
the magnetothermal effect for several special cases, i.e. 
the XY (A = 0) limit, and the low-temperature limits, 
where the analytical solutions arc available. 
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4.1 XY model 

For the XY model, both the spin and thermal currents 
are conserved and hence the magnetothermal effect is 
exactly calculated by using the Jordan- Wigner transfor- 
mation (see Appendix). In Fig. 2, the magnetic Seebeck 
coefficient S together with the thermomagnctic figure of 
merit ZT are depicted for various magnetic fields. In this 
case, one finds the magnetic Seebeck coefficients are al- 
ways negative (S < 0) at any finite temperatures, which 
implies the carriers of the transport are dominated by 
the down spins. 




Fig. 2. Magnetic Seebeck coefficient S (a) and figure of merit ZT 
(b) for the XY model A = 0. For h = h c = 2J, S and ZT 
converge to S = -1.89738 and ZT = 2.02697, respectively (see 
eq. (A-5)). 




Fig. 3. Low-temperature behavior of the magnetic Seebeck coef- 
ficient S in units of T/J for critical regime —1 < A < 1. Inset: 
The magnetic field hi where the low-temperature asymptotics 
changes sign is shown as the function of A. 



to the argument as in rcf. 36: 

S = ~^ forT>J. (4.1) 

From this relation and the results in Fig. 3, we expect 
S has a positive peak and changes sign at certain tem- 
perature T at least for A > 0.5 and h < h\. Namely 
the crossover from the regime dominated by the "down- 
spin" -like carriers to the regime dominated by the "up- 
spin"-like carriers occurs at T . 

In the next subsection, beyond the limitation of the 
analytic methods, we numerically calculate the tempera- 
ture dependence of the magnetic Seebeck coefficients for 
various magnetic fields and anisotropies. 



The magnetic Seebeck coefficient S for h < h c = 2 J is 
linear in T at low temperature. One also observes that 
S diverges at T = as S ~ const. jT (see eq. (A-5) in 
detail) for h > h c , due to the mass-gap (spin insulator). 
Correspondingly the thermomagnctic figure of merit ZT 
increases with increasing the magnetic fields and diverges 
as ZT ~ const./T 2 at T = for h > h c . 

4-2 Low-temperature asymptotics for — 1 < A < 1 

Next we mention the low-temperature behavior T <C 
h < h c in the presence of the interaction strengths — I < 
A < 1. The phenomenological relation yields the leading 
low-temperature magnetic Seebeck coefficient for T <C 
h < h c : 30 

S = a(h)T + 0(T 2 ). 

In Fig. 3, the coefficient a(h) is depicted as a function 
of the magnetic field for various anisotropy parameters. 
For weak interaction strengths A < 0.5, the leading be- 
havior is negative. On the contrary, for A > 0.5, the 
low-temperature behavior changes sign at certain mag- 
netic field h = h\ (sec the inset in Fig. 3). The value h\ 
shifts to higher values with the increase of the interac- 
tion strength. On the other hand, the magnetic Seebeck 
coefficient at high temperatures is determined according 



4-3 Exact diagonalization 

We next present our numerical results obtained by the 
exact diagonalization (ED) method. Using the conserva- 
tion of the total magnetization A4 and the translational 
invariance of the system, we performed a full diagonal- 
ization up to L = 18. To this end, we reduced the integral 
representation of the spin Drude weight Ds in (3.7) (see 
also (2.3) and (2.2)) to a useful form. The spectral de- 
composition H\m) = E m \m) directly yields 

D\ = P 22 Pm\(m\Js\n)\ 2 , Pm = - ^ , 

E m = E n ^ n 

(4.2) 

Applying the partial integration to iss (2-2), one also 
derives 

n" / v\ \( m \Js\n)\ 2 

where K is the kinetic term K = J^2 k (o-^a^ +1 + h.c). 
Here we have used the relation 

/ (Js(-iX)J s )dX = (-K), (4.4) 
Jo 

which is valid in the thermodynamic limit L — > 00. In 
Fig. 4, the size dependence of the magnetic Seebeck co- 
efficient S 1 (S n ) calculated by using D l s (D"). S 1 and 
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14=1.0, h=l ' s i 

5 n 
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*L=\1, 14, 16, 18 



2 4 6 8 10 

T/J 

Fig. 4. Size dependence of the magnetic Sccbcck coefficient S for 
A = 1 and h/J = 1 

S u show a good agreement at high temperatures, con- 
firming Dg (4.2) and (4.3) are equivalent at high 
temperature limit. In this region, one also observes a rel- 
atively weak size dependence. On the other hand, at low 
temperatures, both S 1 and S 11 exhibit strong size de- 
pendences, which mainly stem from the strong finite size 
correction of the spin Drude weight. In this region, one 
sees S l (S u ) decreases (increases) with increasing the 
system size. From this observation, we expect the mag- 
netic Seebeck coefficient S in the thermodynamic limit 
converges to an intermediate value of S l and S u . 

The magnetic Seebeck coefficient S for several 
anisotropics A and magnetic fields h are shown in Fig. 5. 
One observes that S strongly depends on h and A. As 
expected in eq. (4.1), S is negative at high temperatures, 
and converges to zero at T — > oo. 

Above the critical field h > h c where all the spins 
point up at T = 0, the behavior of S is similar to that of 
the XY model. Namely S monotonously decreases with 
increasing h. At low temperature regime, S converges to 
a certain finite value for h = h c , and diverges like S ~ 
—5/T for h > h c , where S is the one-magnon excitation 
gap (see later discussion). 

For h < h c , in contrast, the behavior of the magnetic 
Seebeck coefficients in the regime A > is clearly dif- 
ferent from that for A < 0. Namely for A > 0, S has 
a positive peak (in other words S changes sign) below 
certain magnetic field h < hi, while for A < 0, S is 
always negative as in the XY model. This behavior is 
consistent with the former prediction mentioned in the 
preceding subsection. In our investigation, however, we 
find that this sign changes occur even in < A < 0.5, 
which cannot be concluded as long as we concentrate on 
the low-temperature behavior since the low-temperature 
asymptotics is negative in this region (cf. Fig. 3). One 
can also find in Fig. 5 for A > that the positive peak 
grows in height, and the temperature range where S > 
becomes wide with the increase of A. 

In Fig. 6 the boundary on which S changes sign is 
shown as a line on the h — T plane. The inside (out- 
side) of the curve denotes the region S > (S < 0). 
For any A > 0, the line has a maximum at a certain 




T/J 



Fig. 5. Temperature dependence of the magnetic Seebeck coeffi- 
cients S for both A > and A < 0. 

finite temperature. This feature together with the low- 
temperature asymptotics shown in Fig. 3 indicate that, 
for hi < h < hi, the sign change occurs at the tempera- 
tures Ti and T 2 (S(T) > for T x < T < T 2 and S{T) < 
for < T < Ti or T > T 2 ), while for < h< hi, it oc- 
curs at only one point T = T a (S(T) > for T < T 
and S(T) < for T > T ). Here h 2 is the height of the 
maximum of the line depicted in Fig. 6 and hi is defined 
in the preceding section (see Fig. 3). 

This remarkable feature may be qualitatively inter- 
preted by considering the following extreme cases, (i) For 
fe> 1 where the almost every spin points up, the elemen- 
tary excitation is described by the down-spin magnon. In 
this case, the mobility of the down-spin magnon may be 
larger than that of the up-spin magnon, implying the 
negative magnetic Seebeck coefficient S < 0. (ii) For 
A 3> 1 and < h <C 1, the excitation from the ground 
state is mainly characterized by the up spinon, which also 
plays a role as a carrier, indicating the positive magnetic 
Seebeck coefficient, i.e. S > 0. (iii) For T>1 where the 
interaction A is irrelevant, the behavior of S is similar 
to the XY case where S < 0. 
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4=0.5, 0.75, 1.0, 1.25, 1.5 




S>0 



10 



Fig. 6. Boundary where the magnetic Seebcck coefficient S 
changes sign. 



For the realistic regimes A, T, h ~ 1, the sign changes 
of S might be interpreted as the competition of the above 
three mechanisms. Namely the sign change at T = To or 
T = T2 is explained in terms of (ii) and (iii). On the 
other hand, the sign change at T — T\ is described by (i) 
and (ii), i.e. the crossover from the regime dominated by 
the carriers of down magnons to the regime dominated 
by the carriers of up spinons occurs at T = 7\ . 

We next investigate in detail the divergent behavior of 
S above the saturation field h c . As shown in Fig. 7, the 
divergence of S at low temperature can be canceled by 
multiplying it by T/J. The value of ST/ J approaches 
—5/ J in the low-temperature limit, where 5 is the one- 
magnon excitation gap 5 = h — 2J(A + 1). Thus we 
propose that the leading asymptotics at low temperature 
is written as S — — 5/T. 



"3 




h/J-2(A+l)=1.0, 2.0 



4 6 
T/J 



10 



Fig. 7. The magnetic Seebeck coefficient S multiplied by T/J 
above the critical fields. In the low-temperature limit, S be- 
haves as S ~ — <5/T where <5 is the one-magnon excitation gap 
5 = h- 2J(A + 1). 



For h = h c where S converges to a finite value in the 
low-temperature limit, we observe a universal tempera- 
ture dependence of S for A > as shown in Fig. 8(a), i.e. 
the behavior of S is well described by that of the non- 
interacting case. The existence of such a universal behav- 



ior at h = h c has also be pointed out by Heidrich-Meisner 
et al. in the investigation of the thermal Drude weight. 29 
For A < in Fig. 8(b), in contrast, the temperature 
dependence of 5* strongly depends on the anisotropy. 
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Fig. 8. Temperature dependence of the magnetic Seebcck coeffi- 
cient S for the critical field h c = 2 J(l + A) 



Finally we discuss our numerical results for the ther- 
momagnetic figure of merit ZT. The temperature depen- 
dence of ZT obtained by using Dg 1 is depicted in Fig. 9 
for various magnetic fields and several anisotropics. For 
h < h Cl the figure of merit ZT approaches zero at zero 
temperature. On the other hand, ZT converges to a cer- 
tain finite value for h = h c , and diverges for h > h c , 
which is similarly observed in the XY case. In general, 
for fixed magnetic field, ZT in A > tends to decrease 
with the increase of the interaction strength A, reflecting 
the fact that the thermal Drude weight increases, while 
the spin Drude weight decreases with increasing A. 

5. Summary and Discussion 

In this paper, the magnetothcrmal effect was studied 
for both the ferromagnetic and antiferromagnetic regimes 
of the XXZ chain with external magnetic fields. The mag- 
netic Seebeck coefficients and the figures of merit were 
numerically calculated by using the exact diagonaliza- 
tion up to 18 sites. For the antiferromagnetic regime, we 
found the magnetic Seebeck coefficient changes sign at 
certain temperatures, which is interpreted as an effect of 
the competition among the interaction strengths, tem- 
peratures and external magnetic fields. 

An extension to intcgrablc correlated electron systems 
with the spin degree of freedom such as the Hubbard 
model and the supersymmetric t — J model is an inter- 
esting problem. Due to the integrability, the thermal and 
the spin Drude weights are considered to be finite at finite 
temperatures. In particular, the energy current for the 
supersymmetric t — J model is conserved and hence our 
method developed in this paper can be directly applied. 
For zero magnetic fields, Zemljic and Prclovsek have re- 
cently investigated the thermoelectric power of the Hub- 
bard model. 37 Because of strong correlations, the sign 
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Fig. 9. Temperature dependence of the figure of merit for ZT 
both A > and A < 0. 



changes of the (electric) Seebeck coefficient have been 
observed as in the present case. In the presence of ex- 
ternal magnetic fields, the magnetothermal effect arises 
and crucially affects the charge transport and the ther- 
moelectric effect. The investigation of this effect is also 
of importance for better understanding of the transport 
properties of strongly correlated electron systems with 
reduced dimensions. 

Acknowledgment 

We would like to thank A. Klumper and M. Oshikawa 
for stimulating discussions. This work is supported in 
part by a 21st Century COE Program at Tokyo Tech 
"Nanometer-Scale Quantum Physics" by the Ministry of 
Education, Culture, Sports, Science and Technology KS 
is also supported by Grant-in-Aid for Young Scientists 
(B) No. 17740248. 

Appendix: XY Case 

In this appendix, we provide the exact results of the 
transport coefficients and the thermomagnetic power for 
the XY chain, which corresponds to A = in eq. (3.1). 
In this case, the spin current J7s is a constant of motion, 
i.e. \H,Js\ = 0. Hence from eq. (4.2), one sees that the 
spin Drude weight is simply expressed as D$ — fl{Ji)- 



To calculate the correlation functions (J7|) and (JeJs), 
we introduce the partition function Z as 

, yn , x hiTrexp{-/^ + AiJ s + X 2 Je} n 
lnZ(Ai,A 2 ) = • (A-l) 

Taking the derivative of this with respect to the param- 
eters Ai and A2, and then setting (Ai,A 2 ) = (0,0), one 
obtains the desired quantities: 

D s =(3(Ji)=d 2 Xi lnZ(A 1 ,0)| Ai=Q , 



(Ji)=dl 2 lnZ(0,A 2 )| A2=0 , 

(JeJs) = dx 2 d Xl InZ(Ai,A 2 )| Ai=0)A2=0 . 



(A2) 



For actual evaluation of Z(X\, A 2 ), it is convenient to 
transform the system to the spinless fermion model (3.1 f ) 
(note that we set e = 1). Performing the Fourier trans- 
form, one easily obtains 

TL = 2J^n p cosp— /i^(n p — i), 

p v 

Js = — 2jy^n p sinp, J~ E = — 2J 2 ^ n p sin 2p, (A-3) 
p p 
where p = 2-im/L; n e {-L/2 + 1 , -L/2, ■■■ , L/2 - 
l,L/2}. Substituting them into (A-l) and taking the 
thermodynamic limit L — > 00, we obtain 

lnZ(\ 1 ,\ 2 ) = -^ + ±- J\n(l + e^)dp, 
s(p) = — 2 (3 J cos p — 2Ai Jsinp — 2A 2 J 2 sin2p 



(3h. 
(A-4) 



Thus eq. (A- 2) yields 

27T J_„ 



(3 J 2 sin 2 p 



(Ji) 



1 r 



cosh 2 ((3 J cos p- (3h/ 2) 
J 4 sin 2 2p 



cosh 2 (/3J cos p- (3h/2) 



dp, 
dp, 



(JeJs) 



1 



J 3 ship sin 2p 



-dp. 



2ir J_ v cosh 2 ((3J cosp — 0h/2) 

Combining the above equations with (3.6), (3.9) and 
(3.10), one can calculate the thermal Drude weight and 
the magnetic Seebeck coefficient for arbitrary tempera- 
tures and magnetic fields. The low-temperature asymp- 
totics of the magnetic Seebeck coefficient and the figure 
of merit are explicitly written in the following form. 

7T 2 h 



3v 2 J 
.9(3) 



0(P~ 3 ) for h < h c 



Oi/3- 1 ) forh = h c ' 



5(1) 

(2J- h)P + 0(l) ioih>h c 



3v 2 (3- 



o(r 4 ) 



ZT = I 



for h < h c 
Oifi- 1 ) iorh = h c 
(2 J - h) 2 (3 2 + 0(f3) for h>h c 



5 2 (3) 



.g(l) 5 (5) -.g 2 (3) 
2 
3 



(A5) 
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where v s = \/4</ 2 — h 2 is the velocity of excitations, 
h c = 2 J and g(n) := (1 - 2 1 -"/ 2 )r(l + n/2)((n/2) (note 
that r(x) is the Gamma function and ((x) the Riemann 
zeta function) . On the other hand, the high-temperature 
asymptotics are given by 



l) 

2) 
3) 
4) 
5) 

6) 

7) 
8) 
9) 

10) 

11) 
12) 

13) 
14) 
15) 



-h/3 + 0(f3 2 



ZT~ k2 



(A-6) 
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